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ABSTRACT 

The study of magnetic connectivity in the solar corona reveals a need to generalize the field line 
mapping technique to arbitrary geometry of the boundaries and systems of coordinates. Indeed, 
the global description of the connectivity in the corona requires the use of the photospheric and 
solar wind boundaries. Both are closed surfaces and therefore do not admit a global regular 
system of coordinates. At least two overlapping regular systems of coordinates for each of the 
boundaries are necessary in this case to avoid spherical-pole-like singularities in the coordinates of 
the footpoints. This implies that the basic ch aracteristic of mag netic connectivity — the squashing 
degree or factor Q of elemental flux tubes (jTitov et al.l 120021 ) — must be rewritten in covariant 
form. Such a covariant expression of Q is derived in this work. The derived expression is very 
flexible and highly efficient for describing the global magnetic connectivity in the solar corona. In 
addition, a general expression for a new characteristic Qj^ which defines a squashing of the flux 
tubes in the directions perpendicular to the field lines is determined. This new quantity makes it 
possible to filter out the quasi-separatrix layers whose large values of Q are caused by a projection 
effect at the field lines nearly touching the photosphere. Thus, the value Qj^ provides a much more 
precise description of the volumetric properties of the magnetic field structure. The difference 
between Q and Q± is illustrated by comparing their distributions for two configurations, one of 
which is the Titov-Dcmoulin (1999) model of a twisted magnetic field. 

Subject headings: Sun; coronal mass ejections (SMEs) — Sun: flares — Sun: magnetic flelds 



1. INTRODUCTION 

The structure of magnetic field is often an 
important factor in many energetic processes in 
the solar corona. This especially refers to the 
topological features of magnetic structure such 
as null points, separatrix surfaces, and separa- 
tor field lines. They serve as preferred sites for 
the formation of current sheets and the corre- 
sponding accumulation of the free magnetic e nergy 
( Sweetlll969l:lBaum fc Bratenah]|[l98Q:,Svrovatskii 



vert into other forms: thermal, radiative and ki- 
netic energy of plasma and accelerated particles. 
This process is considered to be a driving mech- 
anism of ma ny energetic phenomena in the solar 
atmosphere (jPriest fc Forbes! I2OOOI : iParkerl Il979l 
19941 ). 

Over the last decade, it also became clear 
that the geometrical analogs of the separatrices 



that the geometrical analogs 01 the separa,trices 
( Longcope fc Strausd [r994allbl: iTitov et all IT999t 



Titov fc Hormdl2002tlTitov et al.H2002f). the so 



■ noiLlr 1^- — llinnnllT t 7<' 1 Jinncl Called quasi-separatrix layers (QSLs , (^Priest fc Demoulin 

198lHLau fc Fmn||l990|:|Longcope fc Cowley|l996t , r—^ W^a^SL .1 IIiqqrJ 
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20011 ). The magnetic reconnection process induced 
in the current sheets at some critical parameters 
allows the accumulated magnetic energy to con- 



I995I: iDemoulin et al1ll996al Il997l )). have similar 
properties. There are indications that the QSLs 
are proba bly more ubiquito us than the true sep- 



aratrices ( Titov et al. 2002f ). This increases the 
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significance of the problem of determining QSLs 
in a given magnetic configuration. In comparison 
witli the separatrices, the determining of QSLs 
requires a more sophisticated technique, which 
is based on a point-wise analysis of the mag- 
netic field line connectivity. The basic quantity 
in this technique is the squashing degree or fac- 
tor Q of elemental magnetic flux tubes. This 
quantity has p reviously been def i ned for the pla- 



nar geometry (iTitov et al.l Il999t iTitov fc Hornij 
2002l: lTitov et al.ll2002l) . which provides a good ap- 



proximation for describing magnetic structures in 
active regions with the characteristic size smaller 
than the solar radius Rq. 

Such an approximation, however, is hardly ap- 
plicable for a global description of magnetic con- 
nectivity in the solar corona including the open 
magnetic field of the coronal holes. The corre- 
sponding large-scale structure of magnetic fields 
is also of substantial interest for solar physics, es- 
pecially, for understanding solar eruptions. So the 
respective generalization of the above technique 
must allow us to determine Q for the coronal vol- 
ume bounded by the photospheric and solar-wind 
surfaces. This immediately raises technical prob- 
lems, which do not exist in the case of the pla- 
nar geometry. First, both these boundary surfaces 
are closed, and therefore, none of them admits a 
global regular system of coordinates. To avoid a 
coordinate singularity of a spherical-pole type, at 
least two overlapping coordinate systems (coordi- 
nate charts) must be used in this case for describ- 
ing the locations of the field line footpoints on each 
of the boundries. Second, the solar-wind bound- 
ary surface generally cannot be a sphere, but some 
other curvilinear surface whos e geometry depend s 
on the coronal magnetic field ( Levine et al. 1982f ). 
These two requirements of the technique can be 
satisfied only by using a covariant approach to 
the description of Q with the coordinate systems 
that are generally different for each of the bound- 
aries. The derivation of such a covariant expres- 
sion for the squashing factor is one of the goals of 
the present work. 

The second goal of the work is to make an es- 
sential refinement of the squashing factor itself. 
The problem is that the large values of Q may 
be caused not only by squashing of elemental flux 
tubes in the volume but also by a projection effect 
at the boundary surfaces. The latter occurs at the 



field lines which are nearly touching the boundary 
at least at one of the footpoints. This effect, in 
particular, tak es place in the vic inity of the bald 
patches (BPs, (|Titov et al.lll993l) ). which are the 
segments of the photospheric polarity inversion 
line (PIL), where the coronal field lines touch the 
photosphere. When analyzing magnetic connec- 
tivity, it is important to discriminate between the 
projection effect and volumetric squashing. For 
this purpose, we derive a covariant expression for 
the perpendicular squashing factor Q_l, which de- 
scribes the squashing of elemental flux tubes only 
in the directions orthogonal to the fleld lines. 

Sections [2] and [3] present the derivations of Q 
and Q±^ and demonstrate on the examples of pla- 
nar and spherical geometry how to apply these 
general expressions. The difference between Q 
and Q ±^ is considered in detail in S|4lby calculating 
and comparing these quantities for two particular 
magnetic configurations. The obtained results are 
summarized in fj5] 

2. COVARIANT FORM OF THE 
SQUASHING FACTOR 

Consider a plasma-magnetic configuration in 
a finite volume with a smooth boundary of an 
arbitrary shape. It may generally consist of 
two or even more surfaces, — for example, the 
photosphere and the solar-wind surface form a 
boundary for the entire solar corona. Each of 
the two footpoints of a given magnetic field line 
may belong in general to any of these surfaces. 
We will use the designations "launch" and "tar- 
get" for the footpoints and parts of the bound- 
ary surfaces at which the field lines start and 
end up. Let (it^,u^) and [w^ ,w'^) be the sys- 
tems of curvilinear coordinates at the launch and 
target boundaries, respectively. The magnetic 
field lines connecting these boundaries define a 
mapping — > {w^,'uP') determined by some 
vector-function {W^{v} ,v?),W'^{u^ .v?')). The lo- 
cal properties of this mapping are described by 
the Jacobian matrix 



D = 



(1) 



For each field line, this matrix determines a lin- 
ear mapping from the tangent plane at the launch 
footpoint to the tangent plane at the target foot- 
point, so that a circle in the first plane is mapped 
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into an ellipse in the second plane (Fig. [T^). 
The aspect ratio of such an ellipse defines the de- 
gree of a local squashing of elemental flux tubes, 
which means that any infinitesimal circle centered 
at a given launch point is mapped along the field 
lines into an infinitesimal ellipse with this as- 
pect ratio at the target footpoint. This gener- 
alizes a coordinate-free definition of the squash- 
ing factor to the case of curvilinear boundaries, 
whose tangent planes are generally not the same, 
as is in the case of p lane boundaries consid ered 
Titov et all (Il999l) : iTitov fc Horlhgl (|2002l ) and 



Titov et al.l (|2002l ). 

To derive an analytical expression for the as- 
pect ratio of the above ellipse, let us introduce 
first a vector function R{u^,u'^) that describes in 
a three-dimensional (3D) Cartesian system of co- 
ordinates the locations of the footpoints at the 
launch boundary. Then the vectors 



dR 



1,2, 



(2) 



determine at this boundary the covariant vector 
basis tangent to the w-coordinate lines. Thus, 



gik = ei-Sk, l,k = 1, 2, 



(3) 



is the corresponding covariant metric tensor, 
which determines local lengths and angles at the 
launch boundary. The dot here stands for the 
usual scalar product in 3D Euclidean space. 

Using ([2]) and the standard Gramm-Schmidt 
procedure, one can construct an orthonormal basis 



ei 
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\/5TT' 

.912 

V.9511 



ei 



£2- 



(4) 
(5) 



Hereafter g = det [gik] is a determinant of the co- 
variant metric tensor. Now, any point of a circle 
of unit radius in the plane tangent to the launch 
boundary is represented by the vector 



o — cos d ei + sin 62 , 



(6) 



whose angle parameter -d e [0, 2tt). 

Suppose that the vector-function R{w^, w'^) de- 
fines the points at the target boundary, then 



= 



(7) 



(a) 




(b) 




Fig. 1. — Linearized field line mapping of a cir- 
cle into an ellipse between tangent planes at the 
launch and target boundaries (a) and between the 
planes perpendicular to the field line at its foot- 
points (6). In general, different and arbitrary coor- 
dinates (■u^,M^) and [w^ ,w'^) with covariant bases 
(61,62) and (61,62), respectively, are assumed at 
the launch and target boundaries. 
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is the field-line mapping image of o at the tan- 
gent plane of the target boundary, where the cor- 
responding covariant basis vectors, parallel to the 
w-coordinate lines, are 



dR 

dw'- 



(8) 



Hereafter a summation over repeating indices with 
their values running from 1 to 2 is assumed. 

With varying the vector 6 traces in this plane 
an ellipse such that 



= g*,oW = -g* 



2"'^ du'' du^ 
X [e\e\ + 6362 -I- cos2i9 {e\e\ - e^e^) 

+ sm2'd{e\el + e{el)], (9) 

where the asterisk indicates that g*j{u^,u'^) is a 
result of evaluating gij {w^ , w"^) at the target foot- 
point {w\w^) = iW\u\u''),W^iu\u^)). 

After some simple trigonometry and lengthy al- 
gebra using equations (HI)-®, equation ^ is re- 
duced to 



52 = -(iV2 + VAr4„4^2 sin 2^), 
where iV^ and A are determined by 



-g 



Ik 



A = 



g d(u^,v?) 



(10) 

(11) 
(12) 



in which g and g* denote the determinants of co- 
variant metric tensors at the launch and target 
footpoints, respectively. The components of the 
contravariant metric tensor g'*^ can be viewed here 
as elements of the inverted matrix [gik] ^ of the 
covariant metric. The value iS is simply plus an 
additional value, which is independent of and 
whose expression does not matter for the present 
consideration. 

What actually matters is that sin 2'd runs values 
from —1 to -1-1 when o{'d) and d{'ff) are tracing, 
respectively, the above circle and ellipse. The min- 
imum — 1 and maximum -1-1 correspond here to the 
minor and major axes of the ellipse, respectively, 
so that its aspect ratio is 



VN^ - 4A2 



1/2 




(13) 



The large values of this ratio do not differ much 
from its asymptotic value 



g = 7vV|A|. 



(14) 



> 



Note also that Q _ 
tion ([13]) yields Q = Omax 



2, since inverting equa- 

/ Omin + Oniin / Omax and 



Omax/omin > 1- Therefore equation will be 
used as a covariant definition of the squashing fac- 
tor. 

It is evident from the derivation of Q that this 
value is invariant to the direction of field line map- 
ping. Indeed, the inverse mapping implies locally 
that 1/omin is a maximum stretching coefficient 
and l/6max is a minimum shrinking coefficient. 
Such coefficients will coincide with the lengths of 
the major and minor axes of the ellipse obtained 
from a circle of a unit radius due to this inverse 
mapping. Thus, although this new ellipse has dif- 
ferent lengths of axes, their ratio is the same as 
for the previous one, which proves the statement. 
A formal proof of the statement is also not diffi- 
cult to obtain by using the derived expressions of 
iV^ and A i n a similar way as in the case of plane 
boundaries ( Titov et al. The invariancy of 

Q to the direction of field line mapping justifies 
its status of a correct measure for the magnetic 
connectivity. 

Note also that A for a given infinitesimal fiux 
tube is a ratio of its cross section areas at the tar- 
get and launch points. Therefore, since the mag- 
netic flux is conserved along the tubes, this value 
coincides with the corresponding inverse ratio of 
the normal field components, so that 



A - bjb:, 



(15) 



4A2 



where i?„ and B* are normal components of the 
magnetic field to the boundaries at the conjugate 
launch and target footpoints. In practice, the nu- 
merical calculation of A through this ratio is more 
precise than that given by equation (|12p and there- 
fore it should be used for computing A in equation 

m- 

The above mathematical construction is re- 
lated to the Cauch y-Green deformation tensor 
(jMarsden et al.ll2002l ) .mown in the theory of elas- 
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ticity. It can be written in our notations as 



Cki 



(16) 



where {W^{u^ ,u'^),W'^{u^ .u^)) and represent, 
respectively, a finite deformation and covariant 
metric tensor of an elastic two-dimensional body. 
The contraction of the Cauchy-Green tensor with 
a pair of orthonormal vectors e,„ and e„ yields the 
tensor 



(17) 



such that its eigenvalues coincide with the squared 
semiaxes o^^^ and o^^^ of the above ellipse. The 
square root of their ratio defines in accordance 
with (fT5)) and the squashing factor Q. 

It should be emphasized that this analogy is 
possible only in our general approach, where two 
independent systems of coordinates are used for 
describing the location of the conjugate footpoints. 
This allows us to apply coordinate transforma- 
tions only at the launch boundary, while keeping 
the coordinates at the target boundary unchanged. 
With respect to these transformations, the object 
defined by (fTH]) does behave as a covariant second- 
rank tensor. The latter is not valid, however, if 
one global 3D system of coordinates is used for de- 
scribing the entire field configuration and so both 
boundaries are subject then to coordinate trans- 
formations. 

This has only a methodological meaning and 
does not exclude, of course, an application of the 
derived expressions to such particular cases. For 
example, consider a closed magnetic configura- 
tion in the half space > with the global 
Cartesian coordinates {x^,x^,x'^) = {u^,u'^,x^) = 
{w^ ,w'^ , x^) and the photospheric boundary plane 
x^ = 0. The field line mapping is then given by 
, X'^ {x^ , x'^)) . There are no more dif- 
ferences between upper and low indices and con- 
travariant (7'*^ and covariant g*j metrics; the latter 
simply turn into Kronecker symbols 5*^' and 6ij. 
So equations (fTTj) . (fT2|) and p5)) are reduced to 



A = 



dX' dX^ 
dx'^ dx^ ' 
d{X\X^) _ B3 

Br 



9(.t1, a;2) 

as required in this case (|Titov et al.ll2002f) . 



(18) 
(19) 



Consider now a more complicated class of con- 
figurations, where both open and closed mag- 
netic field lines are present. Let the configuration 
be described in one global system of coordinates 
{r,6,(j}), where r — Rq corresponds to the pho- 
tospheric launch boundary, while r — repre- 
sents the target boundary. For the open field lines 
reaching the spherical solar-wind boundary of ra- 
dius i?sw, we put i?* = i?sw, while for the closed 
ones we take i?, — Rq. Thus, — — (j), 
11^ = 110^ = and the field line mapping is 
($(0, 9), <d{(t>, 0)), which yields 



and 



sm' 




e 



;^sin ^( 











Rr. 



(20) 



(21) 



where the contravariant metric g'*^ at the launch 
boundary is obtained from the corresponding co- 
variant metric by inverting it simply as a 2 x 2 ma- 
trix. Using these expressions and equation (jlip . 
we obtain 



4 



sin 6 d<^ 

sin 6* d(t) 



1 ae 



sm fc) — — 



99 



(22) 



Equations (|12p and in this case become 
sine 9($,e) Br 



A 



i?| sin( 



d{<p,e) 



— ■ (23) 
B* 



The obtained expressions for N"^ and A have 
seeming singularities at the poles, where they ac- 
tually reduce in the generic case to resolved inde- 
terminacies with $(0,0) and 0(0,0) proportional 
to sin0. These indeterminacies are artificial and 
unrelated to some special properties of the mag- 
netic structure; they are caused by the pole singu- 
larities inherent in the used global spherical sys- 
tem of coordinates. Moreover, their appearance 
is unavoidable in any other global system of co- 
ordinates on a closed sphere-like surface because 
of its intrinsic topological properties. Therefore 
this may generally reduce precision of a numerical 
evaluation of Q near the pole- like points. To avoid 
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such indeterminacies, at least two overlapping co- 
ordinate charts on each of the spherical boundaries 
are required. For this purpose, it is sufficient to 
use two systems of spherical coordinates turned 
with respect to each other on the right angle in 
the ^-direction. Switching from one of such sys- 
tems in its polar regions to the other , as su ggested, 
for example, bv iKagevama fc Satol ( 2004 ). makes 
it possible to resolve the problem of the pole in- 
determinacies. The required expressions for such 
calculations of Q can be obtained again from equa- 
tions (fTT|) and p2|l with properly modified metric 
tensors. 

3. PERPENDICULAR COVARIANT 
SQUASHING FACTOR 

To find the perpendicular squashing factor Qj^, 
we need to know the field line mapping between 
infinitesimal planes orthogonal to the field lines 
at the conjugate footpoints (Fig. Wp)- Note first 
that the projection effect at the boundaries is lo- 
cal, because it depends only on the orientations of 
the tangent planes at the boundaries with respect 
to the vectors of magnetic field at the footpoints. 
So the required mapping between the indicated 
orthogonal planes can be obtained from the re- 
spective mapping between the tangent planes by 
correcting it only at such footpoints. This implies 
that can be expressed in terms of the same 
values as Q and, in addition, the field vectors B 
and B, respectively, at the launch and target foot- 
points. 

We will derive Q±^ by using the same procedure 
as for Q while modifying it in accordance with the 
above comments. The vector o tracing the circle 
of unit radius is given by the same expression 
However, since it lies now in the plane perpen- 
dicular to the field line at the launch point, the 
corresponding orthonormal basis is chosen to be 
orthogonal to B, so that 



ei 



62 



Bxei 

Bxei 
l-BxeiT 



(24) 
(25) 



This vector o is mapped along the field lines into a 
vector o lying in the plane perpendicular to the lo- 
cal field B at the target footpoint. The respective 
mapping dW^ can be represented as a composi- 



tion P ^ o dW o P of three others according to the 
following diagram: 



o 
o 



o 

u 

6 



(26) 



Here the mapping P projects the vector o along B 
onto the plane tangential to the launch boundary 
to yield the vector 

o = o-j^^l^n (27) 

which has a vanishing component along the vec- 
tor eiXe2 perpendicular to such a plane. Then 
this vector o is mapped by the differential of the 
field line mapping dW determined by the Jacobian 
matrix ([1]) into the vector 



(28) 



which lies in the plane tangential to the target 
boundary. Finally, the obtained vector o is pro- 
jected by P~^ along B at the target footpoint onto 
the plane perpendicular to B to result in 



B2 



(29) 



Eliminating now 6 and o from ([27l) ~ ([29)) and 
using ([U with the basis from equations ([24|l -([25 |l . 
we express o in terms of e-basis and magnetic field 
vectors at the conjugate footpoints. This allows us 
to calculate 6^ and then Q± in a similar way as 
done before when deriving Q. The result is 



Q± = 
nI = 

A, = 



Nl/\A, 



f^ d{W\W^) 



Uk 



(30) 
(31) 

(32) 



where the asterisk has the same meaning as in 
equations dS]), HI]) and Thus, the obtained 

Q±_ differs from Q only in the form of the metric 
tensors, which are determined now by 



9Li] - 

Uk _ 



9i3 



9 



Jk 



_ BjBj 
(B-ei X e^f 



(33) 
(34) 
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The asterisk here imphes automaticahy that the 
corresponding values refer to the target footpoint, 
therefore the tilde used for indicating this fact in 
intermediate equation (|29p is omitted in the final 
expression ([55)1 . One can also check that [5"'""^] = 



[5. 



±lk\ 



and 



dot [g 



B2 



so that equation ([32]) reduces to 



|Ai.| 



IB* 



(35) 



(36) 



if, in addition, the magnetic flux conservation is 
taken into account. This expression should be 
used instead of ([32)1 for computing |A_l| in ((30l) 
for the same reason that (|15p should be used in 
([M]). Thus, formulas (I30l)-(l36l) completely define 
the covariant expression for Q^. 

Let us see now how these formulas work in 
the case of a closed magnetic configuration de- 
scribed in a global Cartesian system of coordinates 
{x^,x^,x^) = (u^u^jX^) = (w^,v?,x^) with the 
photospheric boundary plane a;'^ = 0. The formu- 
las are significantly simplified in this case to yield 




(37) 



where the values of the corresponding covariant 
and contravariant components of vectors and ten- 
sors do not differ from each other. 

Note also that expression l|37|) apparently di- 
verges near the PIL, where the normal field com- 
ponent -83 vanishes. This is actually not a true 
singularity but rather an indeterminacy, which is 
resolved to give a low limit of Qj-|pjL = 2 if the 
PIL has no BPs. At the BPs, such an indeter- 
minacy is also resolved but it may generally have 
different limits of Q\_ > 2 at the left and right sides 
of BPs. Thus, Q ^ may experience a jump when 
crossing BPs, unless the configuration is symmet- 
ric as in the example of H4.2I 

Similar to fJH consider also a more general class 
of configurations, where both open and closed field 
lines arc present and bounded by a spherical solar- 
wind surface of radius i?sw and the photosphere 



of radius Rq as before. It is convenient in this case 
to use matrix notations, in which equation (|3ip is 
written as 



where 



D 



/ 5$ a$ \ 

de de 

yikj^ ~d0 J 



(38) 



(39) 



is the Jacobian matrix of the field line mapping. 
The covariant and contravariant metrics at the 
target and launch boundaries, respectively, are de- 
termined by the following matrices: 




where i?, = Rq if a given footpoint belongs to a 
closed field line and i?* = i?sw otherwise. 

The matrix G"*" has two types of singularities, 
which actually lead only to indeterminacies in ex- 
pression (l38|) . The first indeterminacies take place 
at the poles of spherical coordinates {6 = 0, tt). 
They are already discussed above in connection 
with equation ((22| . The second indeterminacies 
take place at the PIL, where Br = . They are 
resolved in much the same way as occurred in the 
previous case of Cartesian geometry to give a fi- 
nite value of Q±|pjj^. If the PIL has no BPs, the 
length of the field lines vanishes near the PIL, so 
that D ^ I and G*j_ G""-"^, which results in 
N'j_ 2. The presence of BPs implies, however, 
a strong projection effect in their neighborhood 
and the corresponding singularities in the Jaco- 
bian matrix D at the BPs. These singularities are 
exactly of the same type and value as those in G^ 
but opposite in signs. So they cancel each other 



in equation p8|) to generally provide different lim- 
its at the different sides of the BPs, as discussed 
above for the case of the plane boundaries. 

It should be also noted that there is one more 
type of plausible singularities not yet discussed. 
These singularities have to appear if a given con- 
figuration has null points of magnetic field and the 
corresponding separatrix field lines. In this case, 
the elemental flux tubes enclosing the separatrix 
field lines split at each of the null points to pro- 
duce singularities in the derivatives of the Jaco- 
bian matrix. Such singularities are due to volu- 
metric rather than surface properties of magnetic 
configurations. Therefore, if present, they appear 
in both Q and Q± distributions thereby indicat- 
ing the existence of magnetic nulls in the corona. 
Since the numerical derivatives are estimated as a 
ratio of finite coordinate differences, their absolute 
values may not exceed the ratio of the coordinate 
range at the target boundary to a chosen incre- 
ment of coordinates at the launch boundary. The 
presence of this upper bound on possible values of 
the numerical derivatives prevents an overflow er- 
ror in computations of the squashing factors at the 
footpoints of the null-point separatrices. A more 
detailed consideration of this type of singularity 
requires a special study, which goes far beyond 
the scope of the present work. 



the other. Such a property is typical for hyper- 
bolic flows in the theory of dynamical systems 
( Arnold IjjS S'l . wh ich was coined in the term HFT 
iTitov et alTf2002l) . The examples below demon- 
strate that Q± is a more accurate quantity than 
Q for characterizing HFTs, although Q also pro- 
vides valuable complementary information on the 
magnetic structure. 

4.1. The simplest possible hyperbolic flux 
tube 

The most simple magnetic configuration in 
which one would expect the presence of an HFT 
is the so-called X-line configuration, whose field 
B is determined in Cartesian coordinates {x,y, z) 

by 



B = {-hx,hy, 1) = V^xVry, 



(43) 



where h = hL/B\\ is a dimensionless field gradi- 
ent characterizing the strength of the transverse 
field JiL compared to the longitudinal field on 
a characteristic length scale L. The right-hand 
side of (I43|l represents the same field in terms of 
the Euler potentials 



V 



ye 



-hz 



(44) 
(45) 



4. COMPARISON OF Q and 

To see the difference between Q and their 
distributions are compared below for several mag- 
netic configurations. An emphasis is made on their 
potentiality for determining so-called hyperbolic 
flux tubes (HFTs). They are deflned as two in- 
tersecting QSLs with extended and narrow pho- 
tospheric footprints characterized by very large 
values of Q. Neglecting the possible curving and 
twisting of an HFT, its cross section variation in 
the longitudinal direction can be respresented as 



(42) 



In other words, the width of one of the QSLs 
is shrunken to the thickness of the other in the 
process of the mapping of their cross sections 
along magnetic fleld lines. This is possible be- 
cause the field lines in HFTs exponentially con- 
verge in one transversal direction and diverge in 



which are constant along the field lines. The use 
of ^ and 77 significantly simplifies the calculation of 
the squashing factors in this configuration. Sup- 
pose that its volume is restricted by \z\ < 1, so 
that the planes z — ±1 are the corresponding 
boundaries. The constancy of ^ and 77 means that 
the boundary points (x_,?/_) and (x+,?/+) are re- 
lated by x-e"'' = x+e'' and y-e'^ = y+e~'^. In 
terms of the notations used in expressions (fT8)) for 
the plane boundaries, this means that Xi = e^^xi 
and X2 = e~'^^X2- Since = 1 everywhere, 
equation reduces simply to |A| = 1, so that 
the resulting squashing factor is 



Q = 2cosh(4/i). 



(46) 



Thus, Q is constant over the entire planes z = ±1 
and, hence, it does not determine any QSL in the 
X-line configuratiq r i with the plane boundaries. 
Priest fc DemoulinI (|l995l ) has arrived at the same 
conclusion by using the norm N defined by (|18p . 
since i?^ = 1, |A| = 1 and so = ^/Q in this 
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configuration. However, they have found with the 
help of N some evidences of the presence of QSLs 
in the X-hne configuration, if it is bounded by cu- 
bic, hemispheroidal and spherical surfaces. 

On the other hand, it is clear from the gen- 
eral point of view that there should be an HFT 
in such a configuration irrespective of the shape 
of the boundary surfaces. Indeed, this configu- 
ration has in the unbounded space two genuine 
separatrix planes x = and y = 0, which can be 
regarded as a limiting case of the bounded config- 
uration with the boundaries moved off to infinity. 
Taking into account also that ft, ^ i, one could 
expect that the proper measure for QSLs must be 
growing with h near the planes a; = and y = 
much stronger than in the remaining volume to in- 
dicate the corresponding QSLs near these planes. 
This would provide in the limit of large h an ex- 
pected continuous transition of such bounded con- 
figurations with QSLs to the unbounded config- 
uration with the genuine separatriccs. The fail- 
ure of Q in determining an HFT in this simple 
case looks very surprising in light of its remark- 
able success in other more complicated field con- 
figurations (see m.2\i . The reason for this failure 
lies actually in the above-mentioned projection ef- 
fect, which is extremely large for the chosen type 
of boundaries. The transverse component here 
grows linearly with the distance from the X-line 
X = y = 0, while the longitudinal component re- 
mains constant. So the farther a given field line 
meets the boundary from the X-line, the more they 
become aligned with each other. 

Such an explanation is fully confirmed by calcu- 
lations of the perpendicular squashing factor Q± 
in this case. To derive an analytical expression of 
Q±, let us choose (^,77) as coordinates on both 
boundary surfaces, so that = = and 
V? — w"^ — -q in (l30 |) -([36 l) . Assume for generality 
that the boundaries are defined by z = ^^±(^,77), 
then according to ([33]) and the vector func- 
tions 

m^r^) = {^e-'^''-,r^e'^^-,Z^), (47) 
- (Ce-"^+,?7e"^+,Z+) (48) 

define, respectively, the launch and target bound- 
ary surfaces. These formulas are needed for calcu- 
lating 0, dSl), dH), §^ dSSl), jHl, dSSl) and (gSl) 
in the chosen coordinates {^,rj). With the help of 



such calculations, equation ([30)1 yields 

^ 2 cosh [2h{Z+ - Z^)]+ (Xl +Y^+Xl+ y?) 

[1 + ft2 {XI + Y^)] [I + ft2 {XI + Y^)] 

in which 

X± - ee-'*^±, (50) 
Y± = rje''^^ (51) 

determine (x, y) coordinates of the conjugate foot- 
points at the defined boundaries. 

These expressions determine Q±^ as a function 
of the Euler potentials ^ and 77. Note also that 
^ — X and r] = y at z ~ 0, so expressions (j49|)- 
()5ip define in addition the distribution of Q± in 
the plane z = 0, to which Q± is mapped along the 
field lines from the boundaries. More generally, 
the combination of these expressions with equa- 
tions ([33]) and (|^5)) . resolved with respect to x and 
y as 

X = ^e-''^ (52) 
y = ^7e''^ (53) 

provides a parametrical representation of Q±, with 
^ and 77 as parameters, in any plane z — const 
between the boundaries. 

In the particular case of plane boundaries, we 
have to put ^■±(^,77) = ±1 in ([31]) -([H]). Using 
then ([44]), ([45]), ^ and ([51]), expression ([49]) for 
Q± can be rewritten even as an explicit function of 
{x,y,z). By comparing ([^ and ([iS]) one can see 
also that the terms containing X± and Y± are re- 
sponsible in this case for eliminating the projection 
effect. As a result of this, the distribution of Q± 
in the z — plane shows very pronounced "ridges" 
along the x and y axes (Fig. [5^) by revealing an 
expected HFT with a characteristic X-type inter- 
section of QSLs along the X-line. According to 
equations ([5^ and ([55]) . such a structure shrinks 
and expands along these axes exponentially fast 
with z to give a typical HFT variation of its cross 
section (see diagram [42]). With growing h, the 
ridges of the Q± distribution and so the corre- 
sponding QSLs become thinner and thinner by 
extending on larger and larger distances from the 
X-line. Thus, the perpendicular squashing factor 
Q± defines indeed an HFT such that it continu- 
ously transforms in the limit of large h into the 
separatrix planes a; = and y = 0. 





f X 

Fig. 2. — Distributions of logQj. in the plane z = for the X-hne configuration from eq. (|43ll at /i = 2. 
The boundaries are two parallel planes z = ±1 (a) or two hyperbolic paraboloids h{x'^ — y^)/2 — z = ^1 (6). 
In the second case, Q± — Q, since the field lines are strictly orthogonal to the paraboloids. The contours 
in both distributions correspond to Q± = 100. In the second case, such a contour represents also the cross 
section of the HFT by the plane z = 0, so that the field lines passing through this contour form the lateral 
surface of the HFT (c) bounded from the top and bottom by the hyperbolic paraboloids. 
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However, the peripheral field lines in this HFT 
are still nearly parallel to the boundaries z = ±1, 
which causes too strong of a variation of the HFT 
cross section. Aesthetically more pleasing HFT in 
the X-line configuration can be obtained by using 
the boundaries which are orthogonal to the field 
lines. These are iso-surfaces of magnetic potential 
F = h (x^ — y"^) /2 — z having the shape of hy- 
perbolic paraboloids. It is natural to chose them 
passing through the points {x ~ y = 0, z ^ ±1), 
which means that F = =pl for such iso-surfaces. 
This condition using ([50|) and ([STjl yields 

y (^2g-2ftZ±_^2g2/.Z±^_^^^1^0^ (54) 

which is a transcendental equation for the Z± 
functions entering in expression (I49|) . This equa- 
tion is not difficult to solve numerically for given 
5 and r] and to use the respective solution for the 
calculation of Q±. An example of the resulting Q± 
distribution in the plane z = at /i = 2 is shown in 
Figure [Ud. The corresponding HFT with the mag- 
netic surface defined by Q±^ = 100 and boundaries 
F = =f1 is presented in Figure For the chosen 
type of boundaries, Q± — Q, so both squashing 
factors define the same HFT. The magnetic field 
for this HFT has the simplest analytical form; the 
hyperbolic paraboloids are also relatively simple 
boundary surfaces. So we believe that this exam- 
ple provides the simplest possible HFT relevant for 
theoretical studies of basic MHD processes, such 
as magnetic pinching and reconnection, in three 
dimensions. 

4.2. HFT in tvi^isted magnetic configura- 
tion 

The considered X-line configuration with z — 
±1 boundaries is an important example, where 
Q±, in contrast to Q, succeeds in determining an 
expected HFT. However, this example is not rep- 
resentative enough to make a general conclusion 
on the potentialities of Q and Q± for detecting 
QSLs. Because the field lines in such a configura- 
tion behave in a rather artificial way over a major 
part of the boundaries. A better comparison of 
Q and Q± can be done by using a more realistic 
field. 

For this purpose, we have chos en the analytical 
model of a twisted magnetic field (jTitov fc Demoulin 



1999( ). hereafter called the T&D model. It de- 
scribes approximate equilibria of a circular mag- 
netic flux rope, whose interior force-free field is 
continuously embedded into a potential back- 
ground field. The latter is produced by fictitious 
subphotospheric sources consisting of two mag- 
netic monopoles of opposite signs and a line cur- 
rent, all located at the axis of symmetry of the 
rope. The axis itself is placed some depth below 
the photospheric plane and the minor radius of 
the rope a is assumed to be much smaller than 
the major one i?c and the distance between the 
monopoles L. 

To compare Q and Q± in detail, we have com- 
puted their distributions for three sets of parame- 
ters which differ only in values of Rc. Two of these 
values {Rc = 85 and 98 Mm) and all the remaining 
parameters are chosen to be exactly the same as 
in the T&D model. By growing Rc but keeping 
other parameters fixed we imitate an emergence 
of the flux rope from below the photosphere. In 
this process, the configuration passes continuously 
through three distinct topological phases. For suf- 
ficiently snMJI__Riij_therejsa^^ BP separatrix 
surface ( Titov fc Demouliiilll999f) — the configura- 
tion with Rc = 85 Mm represents one of these 
topological states. The corresponding distribu- 
tions of Q and Q± are shown in Figures [3^ and 
[Sji, respectively. With growing Rc, this BP and 
the associated separatrix surface bifurcate into two 
parts to give birth to a BP separator field line 
(|Titov fc Demouhnlll999l ) — the configuration with 
Rc = 98 Mm represents this second topological 
phase. The corresponding distributions of Q and 
Q± are shown in Figures [Sb and[3l3, respectively. 
The points Sa and Sd on these figures are the foot- 
points of the BP separator, while Sb and Sc are its 
photospheric contact points. 

Further growing of Rc leads to a complete dis- 
appearance of the bifurcated BP and the associ- 
ated separatrix surfaces — the configuration with 
Rc = 110 Mm represents this last phase. The cor- 
responding distributions of Q and Q± are shown in 
Figures[3l: and[3f, respectively. The magnetic field 
at these parameters becomes topologically trivial, 
since its field line mapping is continuous every- 
where. The whole structure can be continuously 
transformed to a simple arcade-like configuration 
with the help of a suitable photospheric motion. 
Thus, with growing Rc or emerging of the flux rope 
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Rc=85Mm Rc=98Mm Rc^llOMm 
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2 4 6 8 

logQ 

Fig. 3. — Distributions of logQ ([a], [b] and [c]) vs. logQj^ {[d], [e] and [/]) at different values of the major 
radius Rc of the flux rope and with other parameters of the model fixed. The contours of the photospheric 
normal magnetic field = 0, ±100, ±200 and ±400 G (and additionally = ±300 G for [c] and [/]), the 
BPs (thick solid segments of the PIL), and footprints of the BP separatrices (dotted thin lines on [a], [b], [d\ 
and [e]) are also shown. The footprints Sa and 5*^ and touch points Sb and Sc of the BP separator field line 
are indicated on panels (&) and (e). 
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the topological complexity of the configuration 
first increases and then abruptly decreases. 

On the contrary, the geometrical distortion of 
the field line mapping gradually increases in the 
configuration during this process. Both squashing 
factors Q and Q±^ continue to grow with in nar- 
row strips of the photospheric plane. It is clearly 
seen from the Q and distributions plotted with 
the same grayscaling (Fig. [3]) that Q > Q± ev- 
erywhere and the difference between Q and Q± 
becomes smaller with growing i?c. The most sig- 
nificant difference between them is seen at the first 
two phases near the BPs and the footprints of the 
associated separatrix surfaces. As previously an- 
ticipated, Q always rises in these regions, while 
Q± does not, except near the contact points Sb 
and Sc of the BP separator (Fig. [3^). The value 
Q± does rise there but only to give birth to a part 
of HFT footprints, which are matured eventually 
in the third phase (Fig. [3f). As concerned with the 
indicated rise of Q, it is mainly caused by the pro- 
jection effect: the field lines which are close to the 
BP separatrix surfaces approach the photosphere 
near the BPs at a small angle to the horizontal, 
which strongly distorts the footprints of the corre- 
sponding elemental flux tubes. Comparing Q and 
Q± at Rc = 110 Mm reveals that the latter is valid 
for the central part of the PIL as well. 

The discussed features of the distributions be- 
come more transparent if one saturates the grey 
shading in the plots at the values > 2.5 (Fig. [5]). 
These new plots show that the QSLs based on 
the Q± distribution are characterized by a thin- 
ner and more uniform thickness. Their footprints 
acquire at Rc = llOMm a clear fishhook-like shape 
in each of the photospheric polarities with maxi- 
mums of Q reaching ~ lO'*. The QSLs rooted at 
such "fishhooks" intersect eac h other by c ombin- 
ing themselves into an HFT (jTitov et al.l [2003). 
This structural feature seems to be very robust, 
because it appears even in twisted configurations 
which a re not in force-free or m agnetostatic equi- 
librium ( Demoulin et al.|[l996bf) . One can see from 
Figure[5K that, except for an essential twisting dis- 
tortion, the cross section of such an HFT varies 
exactly according to diagram (|42|) . 

It has yet to be proved, but it seems to be 
quite natural that this HFT is pinched into a ver- 
tical current sheet below the flux rope by its up- 
ward movement when the kink or torus instabil- 



ity is developed in the conflguration at a suffi- 
ciently large twist of the field lin e s in the rope 



(|T6rok et al.l 120041 : iRoussev et aD l2003h . This 



interpretation is very important for understand- 
ing the properties of sig moidal struct ures in fiar- 
ing configurations (Kliem et al.l 120041 ). Figure [5^ 
suggests that the sigmoids in such configurations 
are simply pinching HFTs illuminated by a hot 
plasma material which appears there due to the 
reconnection process in the above-mentioned ver- 
tical current sheet. This seems to be valid at 
least for the third topological phase of the fiux 
rope emergence, while an additional interaction 
with the photosphere must be involved at the first 
and seco nd topolo gical phas es, wh e re the BPs are 
prese nt ( Titov fc Demouhn 19991 : iFan fc Gibson 
20031 ). Panels (d)— (/) in Figures H and H demon- 



strate that the footprints of the BP separatrix 
surfaces follow very close to the HFT footprints 
emerging gradually with growing Rc. This implies 
the corresponding similarity in the shapes of such 
separatrix surfaces and HFTs. So the explanations 
of the sig moids that rely on e ither the presence of 
the B Ps (iTitov fc DemouHnl[l9 99: Fa n fc Gibson 



20031 ) or the HFTs (jKliem et al.ll2004f ) are not al- 



ternative but rather complementary, since they re- 
fer to different phases of the flux rope emergence. 

It should be noted also that both Q and Q± 
distributions (Fig. [4]) contain at the border of the 
fiux rope two less pronounced horseshoe-like fea- 
tures with maximums ~ 10^. The QSL rooted at 
these "horseshoes" has a helical shape (Fig. [5)d) 
with a slightly varying cross section along the field 
lines. So this QSL has a structure qualitatively dif- 
ferent from those two which form the HFT. The 
comparison of Q and Q± shows that the squash- 
ing of the fiux tubes in this helical QSL is only in 
part due to the projection effect. The Q± distri- 
bution demonstrates that the major contribution 
to Q comes from the shearing of the twisted field 
lines in the rope. Thus, both distributions reveal 
a helical QSL which is a part of the inner border 
layer of the flux rope. In this respect, the con- 
sidered example demonstrates that our squashing 
factors help identify flux ropes themselves. It is 
not a problem, of course, to locate the flux rope in 
the TfcD model, where its parameters are known 
from the construction of the model. Yet identify- 
ing flux ropes in more complicated configurations 
obtained, for instance, numerically from magne- 
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Fig. 5. — HFT (a), the helical QSL (c) and their cross sections, for the HFT (6) and for the QSL (rf), in 
the T&D twisted magnetic configuration at i?c = 110 Mm. Other parameters of the model, as well as the 
contours of the photospheric normal magnetic field B^, are the same as in Figs. [3]: and|3f; the white stripes 
on the photosphere represent the HFT and QSL footprints, whose contours correspond to Q = 100. 
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togram data is a real problem. If such configu- 
rations are topologically trivial, like those in the 
third phase of our example, the determination of 
QSLs seems to be the only method for identifying 
flux ropes. 

As shown above, both squashing factors allow 
us in the case of the T&D model to determine 
similar QSLs, except that Q±_ is more advanced 
than Q near BPs, whenever they appear. There- 
fore, wc think that, in general, if the numerical 
grid used for computing Q±_ is fine enough for de- 
tecting possible BP separatrix surfaces by sudden 
spikes in the distribution, the use of only Qj. 
would be sufficient for the structural analysis of 
configurations. In practice, however, the required 
resolution of the grid cannot be always easily fore- 
seen. Also the computational cost for Q is not 
really high, since the same input data as for Q± 
can be used. So it is sensible to compute both 
these distributions at a time and compare them 
in the same way as we did in the considered ex- 
ample. Some redundancy of the information con- 
tained in such distributions is not superfluous but 
useful, especially, in the case of complicated real 
magnetic configurations. Thus, from the practical 
point of view, the value Q should not be consid- 
ered as obsolete but rather as a complementary 
characteristic of magnetic connectivity. 

5. SUMMARY 

We have derived a covariant form of the squash- 
ing factor Q, which enables us to determine quasi- 
separatrix layers (QSLs) in both closed and open 
magnetic configurations with an arbitrary shape of 
boundaries. The corresponding expression for Q 
assumes that the Jacobian matrix of the field line 
mapping and the metric tensors at the footpoints 
are known. The expression admits also that such 
"input data" can be represented with the help of 
two different coordinate systems for determining 
location of the conjugate footpoints on the bound- 
aries. This provides a firm theoretical basis for a 
global description of the field line connectivity in 
the solar corona. 

To eliminate the projection effect at the field 
lines which are nearly touching the boundary, the 
perpendicular squashing factor Q±_ is also derived 
in a similar covariant form. The value Q± defines 
the degree of squashing of elemental magnetic flux 



tubes only in the directions orthogonal to the field 
lines. In the definition of Q_\_, the boundaries en- 
closing the magnetic configuration constrain only 
the length of the flux tubes while not affecting 
their cross sections at the footpoints. For calcu- 
lating (5_L, the vectors of magnetic field at the foot- 
points arc required in addition to the same input 
data as for Q. The use of both covariant squash- 
ing factors is demonstrated by calculating them 
for the boundaries with the planar and spherical 
geometries. Then the properties of Q and Q±_ are 
compared by considering two examples of mag- 
netic configurations. 

The first example is a classical X-line configura- 
tion of potential magnetic field in a plasma volume 
restricted by two boundary surfaces. It is easy to 
show that for the plane boundaries perpendicular 
to the X-line, the value Q is constant. So the Q 
distribution does not allow us to define any QSLs 
in such a configuration. The reason for this failure 
lies in the projection effect, which is very strong 
for the field lines distant from the X-line. We have 
also calculated an analytical expression of Q±_ for 
the same field but with the boundaries of an ar- 
bitrary shape. In the case of the plane bound- 
aries, this new value Q_\_, in contrast to Q, has 
a non-uniform distribution, which does reveal the 
expected two QSLs. These QSLs intersect each 
other by combining themselves into a hyperbolic 
flux tube (HFT). A more elegant HFT is obtained 
for the X-line configuration with the boundaries 
orthogonal to the field lines. 

To make a better comparison of the proper- 
ties of Q and Q^, a second magnetic configura- 
tion more relevant for solar physics is considered. 
The respective field is defined by using the Titov- 
Demoulin (1999) model of a force-free flux rope 
embedded into a potential background field. Con- 
trary to the case of the X-line configuration re- 
stricted by the plane boundaries, both the Q and 
Qj^ distributions reveal QSLs in the twisted con- 
figuration. These distributions are similar every- 
where except near bald patches (BPs) and foot- 
prints of the associated separatrix surfaces, when- 
ever the BPs exist. By definition, the value Q±_ 
is free of the projection effect, so that Q±_ rises 
near BPs only if the corresponding fiux tubes are 
subject to a volumetric squashing. This is not the 
case, of course, for the value Q, which always rises 
in such regions of the photosphere. So, in compar- 
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ison with Q. the value Q±^ shows itself to be again 
a superior characteristic for analysis of magnetic 
connectivity. 

Nevertheless, we have argued that it is more 
practical in general to compute both squashing 
factors for analyzing the structure of a given mag- 
netic configuration. This docs not require addi- 
tional significant effort, while making it easy to 
discriminate between the volumetric squashing of 
elemental flux tubes and the surface projection ef- 
fect at the boundaries. 
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